Distribution functions in non-equilibrium theory of superconductivity 
and Andreev spectroscopy in unconventional superconductors 
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We present a new theoretical formulation of non-equilibrium superconducting phenomena, includ- 
ing singlet and triplet pairing. We start from the general Keldysh-Nambu-Gor'kov Green's functions 
in the quasiciassicai approximation and represent them in terms of 2x2 spin-matrix coherence func- 
tions and distribution functions for particle-type and hole-type excitations. The resulting transport 
equations for the distribution functions may be interpreted as a generalization to the superconduct- 
ing state of Landau's transport equation for the normal Fermi liquid of conduction electrons. The 
equations are well suited for numerical simulations of dynamical phenomena. Using our formula- 
tion we solve an open problem in quasiciassicai theory of superconductivity, the derivation of an 
explicit representation of Zaitsev's nonlinear boundary conditionsu at surfaces and interfaces. These 
boundary conditions include non-equilibrium phenomena and spin singlet and triplet unconventional 
pairing. We eliminate spurious solutions as well as numerical stability problems present in the orig- 
inal formulation. Finally, we formulate the Andreev scattering problem at interfaces in terms of the 
introduced distribution functions and present a theoretical analysis for the study of time reversal 
symmetry breaking states in unconventional superconductors via Andreev spectroscopy experiments 
at N-S interfaces with finite transmission. We include impurity scattering self consistently. 



PACS numbers: 74.20.-z, 74.80.Fp, 74.50.+r 
I. INTRODUCTION 

Conduction electrons in normal metals are generally 
well described by Landau's Fermi liquid theory.0 Ac- 
cording to Landau a system of strongly interacting elec- 
trons may be viewed as an ensemble of quasiparticles 
which can be represented by a classical distribution func- 
tion and obeys, a semiclassical Landau-Boltzmann trans- 
port equations This semiclassical behavior of a quantum 
many-body system is a consequence of Pauli's exclusion 
principle which restricts the momentum space accessible 
to low-energy quasiparticles to a thin shell near the Fermi 
surface. The ratio of the volume of the accessible mo- 
mentum space to the total volume enclosed by the Fermi 
surface is of the order ksT/Ef <C 1, and is the fundamen- 
tal expansion parameter of Fermi liquid theory. Landau's 
Fermi liquid theory is exact in leading order in an asymp- 
totic expansion in ksT '/ 'Ep and other small parameters of 
an electronic Fermi liquid such as 1/kf^T, hto/Ef, l/kf£, 
where Ef, kf, £t, w, and I are Fermi energy, Fermi wave 
vector, thermal coherence length (£t = hvf /2irkBT), 
frequency of time-dependent perturbations, and quasi- 
particle mean free path. Phase space arguments can 
be used to derive Landau's Fermi liquid theory by con- 
verting a formal diagrammatic expansion of many-body 
Green's functions into ap|-|asymptotic expansion in the 
above small parametersJ3~0 Only a few of the resummed 
Feynman self-energy diagrams contribute in leading or- 
ders, and the dynamical equations for Green's functions 
can be transformed into Landau's transport equation for 
quasiparticle distribution functionsH'Da'cl The price one 
has to pay for the simplifications of the quasiparticle the- 
ory is the need to introduce phenomenological parame- 
ters, such as the quasiparticle velocities and quasiparticle 



interactions In the absence of first principles calcula- 
tions these material parameters have to be taken from 
the experiment. 

The development of a generalized Fermi-liquid the- 
ory for the superconducting state staxted with the work 
by GeilikmanaEl and Bardeen et al£3 shortly after the 
BCS-theory of superconductivity was publishedO These 
authors presented a generalization of the semiclassical 
transport equations of the normal state to the supercon- 
ducting state and used them to explain the electronic 
thermal .conductivity of superconductors. Several early 
workalaiU on transport and linear response in super- 
conductors showed that various semiclassical concepts of 
Landau's Fermi liquid theory could be readily general- 
ized to the superconducting state. A novel semiclassical 
approach to superconductors in equilibrium was initiated 
by de GennesEj who formulated the equilibrium theory 
of superconductivity near T c in terms of classical cor- 
relation functions. These developments established the 
usefulness and accuracy of semiclassical concepts for su- 
perconductors but did not provide a complete semiclassi- 
cal theory of superconductivity. Nevertheless, these early 
semiclassical works were the predecessors_of a compre- 
hensive theory developed by Eilenbergerjl3 and Larkin 
and OvchinnikovEH for superconductors in equilibrium. 
This theory was generalized to non-equilibriumjjhenom- 
ena by EliashbergQ and Larkin & Ovchinnikov.li3 We re- 
gard this theory as the proper generalization of Landau's 
Fermi liquid theory to the superconducting state^ and 
call this theory, following Larkin and OvchinnikovE3, the 
quasiciassicai theory of superconductivity. 

The derivation of the quasiciassicai equations starts 
from Gor'kov's formulation of the theory of supercon- 
ductivity in terms of Nambu-Gor'kov matrix Green's 
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functions Jla Typical spatial variations of the order pa- 
rameter occur on a scale of the coherence length, £o = 
Hvf /2-KkBTc, and the typical time scale, to-, is given 
by the inverse gap, to ~ h/A. Both scales are usu- 
ally large in superconducting metals as compared to 
the corresponding atomic scales, ~ and ~ h/Ef. 

Gor'kov's Green's functions contain detailed information 
on atomic scale properties which average out on the su- 
perconducting scales. To derive the quasiclassical equa- 
tions one has to integrate out atomic scale features in 
the Green's functions, but keep all relevant informa- 
tion for superconductivity. The resulting "^-integrated" 
Green's functions (quasiclassical Green's functions) vary 
on the large scales £o and to and are free of the irrelevant 
fine-scale structures. The quasiclassical equations should 
be compared with Andreev's equations^ which he ob- 
tained by factorizingj-ojat. rapidly oscillating terms in Bo- 
goljubov's equationsEirta Andreev's method is equivalent 
to the quasiclassical theory for superconductors with in- 
finitely long-lived quasiparticles, i.e. without impurities, 
electron-phonon coupling or electron-electron scattering. 
Both theories give identical results in these cases. Hence, 
the quasiclassical theory of Eilenberger, Larkin, Ovchin- 
nikov, and Eliashberg may be considered a generalization 
of Andreev's equations to systems with disorder and fi- 
nite lifetimes of quasiparticles. The generalized theory 
covers basically all phenomena of interest in supercon- 
ductivity. 

One distinguishes in the quasiclassical theory of su- 
perconductivity external (classical) and internal (quan- 
tum mechanical) degrees of freedom of a quasiparticle. 
Quasiparticles move along classical trajectories, and are 
reflected by walls or other obstacles and scattered by 
collisions with impurities, phonons or with other quasi- 
particles. This classical dynamics in coordinate space 
should be contrasted with the quantum-dynamics of in- 
ternal degrees of freedom of a quasiparticle which are the 
spin (S=i) and the particle-hole degree of freedom. The 
quantum coherence of particle and hole excitations is an 
essential feature of BCS pairing theory. It is absent in the 
normal state and is the origin of all typical superconduqt-. 
ing phenomena including the opening of an energy pgap, til 
Andreev's retro-reflection,Ej Tomasch oscillations, Ej vor- 
tex bound states J23 etc. nr-|n 

In a standard notation of quasiclassical theoryOO'E£l 
the distribution functions for quasiparticles are 4x4 ma- 
trices which reflects their quantum mechanical structure 
as density matrices in the 4-dimensional Hilbert space of 
internal degrees of freedom (Nambu-Gor'kov space). The 
matrix elements are functions of the classical variables p / 
(Fermi momentum), R (position), e (energy measured 
from the chemical potential), and t (time), which de- 
scribes the classical degrees of freedom. The external 
motion of a quasiparticle and its internal particle-hole 
state are-coupled in a subtle way as first discussed by 
AndreevE3 

The complexity of dynamical phenomena in supercon- 
ductors makes the elimination of atomic scale fine struc- 



ture an important step towards a solution of dynamical 
problems. The dynamical equations of the quasiclassi- 
cal theory, which are free of microscopic fine structures 
from the outset, can be formulated most compactly by 
using the Green's function technique of Kcldysh.Ea It is 
often useful to distinguish two sources of time-dependent 
phenomena. Firstly, time dependences can arise from 
changes in the occupation of quasiparticle states. In the 
normal phase of a Fermi liquid this is indeed the only 
fundamental dynamical process. Here, the quasiparticle 
states are robust and changes of the quasiparticle wave 
function can be neglected. This is no longer the case in 
the superconducting phase. Quasiparticle states in su- 
perconductors are coherent mixtures of particle and hole 
states determined by the superconducting order param- 
eter. Since the order parameter will, in general, change 
in a dynamical process the quasiparticle states will also 
change. Superconducting dynamics is thus governed by 
the coupled dynamics of both the quasiparticle states and 
their occupation. The Keldysh technique is convenient in 
this case since it works with two types of Green's func- 
tions (g R ' A and g K ) and can be used to introduce dynam- 
ical spectral functions describing the time development of 
quasiparticle states and dynamical distribution functions 
describing the time-dependent occupation of the states. 
Dynamical distribution functions in the superconducting 
state were introduced by Larkin & OvchinnikovEa and by 
Shelankov.Ell 

In this paper we present an exact parameterization of 
the quasiclassical Keldysh Green's functions in terms of 
four coherence functions and two distribution functions. 
The coherence functions are generalizations of the Riccati 
amplitudes introduced recentlyc3nJ for superconductors 
in equilibrium, whereas the distribution functions are the 
generalizations of the distribution function of Landau's 
Fermi liquid theory of the normal state. Compared to the 
conventional quasiclassical theory our formulation leads 
to intuitively appealing and explicit boundary conditions 
at surfaces and interfaces, is numerically very stable, and 
allows for a more transparent interpretation of quasiclas- 
sical dynamics in terms of particle-type and hole-type 
excitations. 

The general framework of the quasiclassical theory 
is briefly reviewed in sections II, III, where we also 
introduce our notation. Dynamical equations for the 
coherence functions are derived in section IV together 
with dynamical equations for the distribution functions 
(transport equations). In section V we solve Zaitsev's 
non-linear boundary conditions for quasiclassical Green's 
functions at interfaces, and obtain physically appealing 
boundary conditions for our coherence functions and dis- 
tribution functions. In section VI we present the gen- 
eral linear response equations in terms of the introduced 
functions. Finally, we formulate in section VII the An- 
dreev scattering problem at interfaces between a normal 
metal and an unconventional superconductor using the 
new theoretical formulation and the resulting boundary 
conditions. We present results for the Andreev reflec- 
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tion amplitudes and the regularly reflected amplitudes 
at (110) interfaces and (100) interfaces between a nor- 
mal metal and a layered d-wave suoetconductor. Our 
formulation generalizes earlier workEJEj to include disor- 
der. We propose an anomalous feature in the reflection 
amplitudes for (110) interfaces as a test for time reversal 
symmetry breaking states. This feature, a strong sup- 
pression of the regular reflection for low-energy quasi- 
particles at interfaces with finite transmission, is sen- 
sitive to sign changes in the order parameter, and has 
the same origin as the zero-energy surface bound states. 
Combined with this suppression is an enhancement of 
the excess current due to Andreev reflection for low en- 
ergy quasiparticles. The sensitivity of this phenomenon 
to time reversal symmetry breaking states provides a new 
tool to study the symmetry of the order parameter. We 
study the effect of disorder on both regularly and An- 
dreev reflected currents. For Andreev spectroscopy in 
unconventional superconductors the low-energy behavior 
of regular reflection is the spectral feature most stable 
against disorder. 



II. KELDYSH SPACE STRUCTURE 

The fundamental quantity in non-equilibrium quasi- 
classical theory of superconductivity is.. the. quasiclassi- 
cal Green's function g = <?(p/, R, e, i)BE-atJ It is a 2x2 
Keldysh matrixES of the form 



9 



6 r 



(i) 



where the elements are 4x4-Nambu matrices, which de- 
scribe the two important residual quantum mechanical 
degrees of freedom: the spin degree of freedom and the 
particle- hole degree of freedom. g R — g R (pf, R, e, t) 
is the retarded, g A = g A (pf, R, e, t) the advanced and 
gK _ g K (p^ R ; e , t) the Keldysh Green's function. The 
classical degrees of freedom are described by a motion 
of the quasiparticles along classical trajectories. All tra- 
jectories through a spatial point R are parameterized by 
the Fermi momentum, pf, and their directions coincide 
with the directions of the Fermi velocities, v/(p/). Along 
a given trajectory with fixed pf all quasiparticles travel 
with the same velocity, v/(p/). In general there can be 
several branches of quasiparticles moving with the same 
velocity but having different momenta. Also the direc- 
tions of pf and Vf(pf) are generally different. However, 
for spherical or cylindrical Fermi surfaces p/ and V/(p/) 
differ only by a scaling factor. 

The quasiclassical Green's function is solution of the 
following transport equation along a given trajectory, and 
of the corresponding normalization conditioroE^HLa (the 
(gi -product is noncommutative and is explained in Ap- 
pendix |X|) 



[e - h, g] gj + ihvfVg = 6; 



Here e represents the energy variable, and h combines 
the molecular (or mean) field self-energies, & m f, the im- 
purity and electron-phonon self-energies, a^, and external 
potentials, v ex t 



h = a m f + <7j + v ext = 



h R h h 
6 h J 



(3) 



which have the following explicit Keldysh structure (f, 
denote Pauli matrices in the particle-hole space) 



eh 

6 ef 3 

a R a K \ 

6 a A r Vext 



_ / cr m/ 







(Tmf 

Vext 6 
6 Vext 



(4) 
(5) 



Disorder will be included by following standard aiierag- 
ing procedure for dilute impurity concentrations.^ We 
denote impurity self-energies by C7 imp . In quasiclassical 
approximation [I 3> 1/kt), the impurity self-energy can 
be written in terms of the concentrations, Cj, of impuri- 
ties of type i and the single impurity t-matrices, ti, 



N 



(Timp (P/ ) 



The t-matrices are solutions of the following equations 
(we suppress the variables R, e, t for convenience) 

U(Pf,p' f ) = "i(p/,P/) + 

iV / (u i (p /) p^)®5(P/)®ti(p/,P/)) , (7) 
\ / Pf 

where Ui(pf,p'^) is the scattering potential of an impu- 
rity of type i. The Fermi surface average (• • -) p ^ is ex- 
plained in Appendix The impurity potential is diago- 
nal in Keldysh space. 



III. NAMBU-GOR'KOV SPACE STRUCTURE 

We parameterize the elements of the Nambu matrices 
in the following way 



g R ' A 



9 RA 


f H,A* 




g R A , 











,9 = 



9 K r 

-f K -~9 K 



12 K 



A K 



(8) 
(9) 



9®9 : 



-7T 2 l 



(2) 



Here g R A , f R,A , A K etc. are 2x2 spin matrices. 

The molecular fields are determined by Landau's quasi- 
particle interaction function, A(pf,p'j), leading to a self- 
energy spin matrix, v m f{pF, R, e, i), which is diagonal in 
particle-hole space. In superconductors this interaction 
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must be supplemented by the pairing interaction of quasi- 
particles, V(p/,p^), which lead to an off-diagonal self- 
energy in particle-hole space, A m f(pp, R, e, t). Thus, 

a m f = P m f + A m f. (10) 

The mean- field-. self energies, Eq. (|l0|), are diagonal in 
Keldysh space.El Their matrix structure in Nambu space 
is 

* f A mf \ ( v m} \ f . 

Am/ = { A mf J ' = { 0% m/ ) ■ C n > 

Not all the matrix elements are independent feom each 
other, but are related by symmetry relations^ For in- 
stance, a quantity x and the conjugated quantity x are 
related by, 

z( P/ ,R,e,i)=z(- P/ ,R,-e,i)*. (12) 

The conjugation operator ( " ) defines an important 
transformation of quasiclassical Green's functions and 
self-energies. We will use it extensively in the following. 

IV. COHERENCE FUNCTIONS AND 
DISTRIBUTION FUNCTIONS 

The numerical solution of the transport equations for 
the quasiclassical Green's functions can be simplified con- 
siderably by introducing a special parameterization in 



terms of 2x2 spin matrix coherence functions, r ) R ' A , j R ' A , 
and distribution functions, x K , and i K , which transforms 
the original boundary value problem for g into initial 
value problems for the coherence and distribution spin 
matrices. The normalization condition is in this formu- 
lation eliminated completely. We present here the re- 
sulting equations and refer for their derivation to Ap- 
pendices |b] and |c[ Before doing this we give a short 
physical interpretation for the coherence functions. In 
the absence of particle-hole coherence, like in the equi- 
librium normal state, the functions 7 fl,A , 7 R,A vanish. A 
superconductor, or a normal metal in proximity to a su- 
perconductor, can be described in equilibrium and in the 
clean limit by Andreev's equationsE3 with Andreev am- 
plitudes u and v. Then, the coherence function "f R , for 
example, is given in terms of the u- and u-spin matrices 
(for positive energies) by the solution of the linear system 
J2a u a/3lps = v aS- Thus, the coherence functions are the 
transformation matrices between the particle and hole 
like Andreev amplitudes. In the presence of quasiparti- 
cle damping the Andreev description breaks down, never- 
theless one can define generalized amplitudes u R,A and 
v R,A . In non-equilibrium they are defined by relations 
like u H (»7 R = v R . Note that these generalized ampli- 
tudes are defined by the quasiclassical Green's functions, 
not by wave functions. 

The quasiclassical Green's functions are conveniently 
parameterized by 



g =- 2ni N {-(^ 9 ^-^ 9 r) (i--7 R »^«7 A ) ) ® N ' (14) 

with the 'normalization matrices' 

\ (l-7«- A ^ 7 «' A )- 1 J ■ [ib > 

In ( |l3| ) the factor N R,A may be written on the left- or right-hand side. The transport equations for the 2x2 spin 
matrix functions are 

ih V/ V 7 R ' A + 2ej R - A = 7 R " 4 s A R ' A 8 1 R ' A + (z R A g 7 fl ' A - 7 R ' A g E R ' A ) - A R - A , (16) 

i?i v f Vj R ' A - 2e^ A = 7 fl - 4 g A R ' A g f-* + (± R A g 7 fl,A - 7*' A g - A R,A , (17) 

ihv f Vx K +ihd t x K + (- 7 R ®A R -?, R )®x K + x K ®(- A A ®7 A + S A ) = 

= -7 R ®E K ®7 A + A A "®7 A +7 n ®A K -S K , (18) 

ihv f Vx K -ihd t x K + (-7 H ®A H -E H )®^ + 5 iC ®(-A A ®7 A + E A ) = 

= -7 K ®£ K ®7 A + A A "®7 A + 7 R ®A K -S K . (19) 
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Equations ([[3]) , (|T^) , and ( |l7| ) generalize a useful formu- 
lation of the equilibrium theory in terms of Riccati-type 
transport equationsr^r 9 ! to noa-equilibrium phenomena. 
Equations (|l3|)-([l9|) are new,E3 numerically very stable 
and provide an efficient way to solve non-equilibrium 
problems in superconductors. Equations (|l^)-(|l9|) need 
to be supplemented by initial conditions. They are im- 
posed for 7 R , 7 A , and x K at the beginning of the tra- 
jectory, and for j A , j R , and x K at the end of the trajec- 
tory. For correctly chosen initial conditions the transport 
equations for j R , j A , and x K are stable in positive vy- 
direction, and the transport equations for 7 A , j R , and 
x K are stable in negative v^-direction. In addition to the 
conjugation symmetries, 



7 R > A ( P/ ,R,M)= 7 ^(- P/ ,R,-M)* , (20) 
*(- P/ ,R,-e,*)*, (21) 



x K {p f ,TL,e,t) 



the coherence and distribution functions obey the follow- 
ing symmetries 



7 A ( P/ ,R,e,t) = 7 R (P/,R,e,i) t 
x K (p f ,IL,e,t) = x K {pf,R,e,t) J < 



(22) 
(23) 



Note that the x K (pf, R, e, t), x K (pf, R, e, t) are her- 
mitean spin matrices. In equilibrium, 



(l- 7 R 7 A )tanh— , 
-(l-fV)tanh^. 



(24) 
(25) 



V. EXPLICIT SOLUTION OF ZAITSEV'S 
NONLINEAR BOUNDARY CONDITIONS 

In the previous sections we have introduced a 
parametrization of the non-equilibrium Keldysh Green's 
function, g, in terms of four coherence functions and two 
distribution functions (2x2 spin matrices) 



9 = gh R ,r,i A ,r,x K ,x K ] 



(26) 



An important problem is the formulation of bound- 
ary conditions— .for these parameters at surfaces and 
interfaces.EieZrEil A boundary condition for g was ob- 
tained by ZaitsevJH which in principle solves this prob- 
lem for perfect interfaces. However, Zaitsev's non-linear 
boundary conditions have unphysical spurious solutions 
which require special care, e.g. in a numerical imple- 
mentation. A linearization of Zaitsev's boundary condj-. 
tions for the equilibirium was achieved recently by YipEil 
for the case of an interface connected to infinite half 
spaces. Our solution generalizes these results to any in- 
terface geometry and to non-equilibrium phenomena. Za- 
itsev's condition relates the quasiclassical Green's func- 
tions with Fermi velocity pointing in direction towards 



r*r, A X K 




1 >"i 

1 1 1 > A 1 



FIG. 1. Notation for the Green's functions at the inter- 
face. Indices 1 and 2 refer to the sides of the interface. The 
arrows for the Fermi momenta (dotted) are for particle like 
excitations. The Fermi velocity directions (full lines) are 
given by the directions perpendicular to the Fermi surface 
(full curves) at the corresponding Fermi momentum. The 
components of the Fermi momenta parallel to the surface 
are conserved (indicated by the thin dotted lines). 

the surface, <7i,j n , g<i in, and those with Fermi velocity 
pointing away, <?i )0 ut, 32, out- Indices 1 and 2 refer to the 
two sides of the interface (see Fig. [j]). Using the defini- 



tions 
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which fulfill the relations 










P a ®P 1 +P 1 ®P a =Q , 


Pa 


%Pa + 


Pi® Pi 


= 1 , 




Pa 


%Pa + 


P 2 ®P2 


= 1, 



(28) 



(29) 
(30) 

Zaitsev's boundary conditions readEl 

[(I -Pa) ®Pl®P2- (I + Pa) ®P2®Px] (1 - H) = 

= -2P ®(l-P ®P a )(l+ft). (31) 

Here, and in the following 1Z = TZ(pj) and V — T>(pj) 
are the reflection and transmission coefficients of the in- 
terface for quasiparticles in the normal state correspond- 
ingly, TZ(pf) + T>(pf) = 1. In the following we present 
the explicit solution of Zaitsev's nonlinear boundary con- 
ditions at spin-conserving interfaces in terms of the co- 
herence functions and distribution functions introduced 
above. It is useful for this purpose to introduce a no- 
tation which indicates the stability properties of solu- 
tions of the transport equations. We use capital letters 
(T R - A , T R - A , X K , X K ) for functions, which are stable so- 
lutions when integrating the transport equation towards 
the surface. Small case letters {j R,A , j R - A , x K , x K ) are 
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used for functions, which are stable in the direction away 
from the surface. We also generalize the notation for the 
conjugation operation. It includes an conversion from 
small case to capital case letters, 

s(p / ,R,e l i) = X(-p / ,R,-e > *)*. (32) 

By integrating in direction towards the surface, the quan- 
tities Jj' A , lf' A i x f:Xf, (j = 1, 2) are known. The quan- 
tities Tf A ,tf A ,Xf,Xf are to be determined by inte- 
grating in direction away from the surface. At the sur- 
face the second set of quantities is determined in terms 
of the first one by boundary conditions. 

The incoming quasiclassical retarded Green's functions 
(with velocity direction towards the interface) on each 
side of the interface are given then by (see Fig. [[]), 

ff Mn = ff[7f,ff,rf,7f,a ; f ! Xf] , (33) 

mn = g[l%, f * , T A , l£,x%,X?} , (34) 



and the outgoing ones (with velocity direction away from 
the interface) 

h,out = ,g[rf,7f ,it ,fi ,x? ,if] , (35) 
h,out = g[ri tT, it f A , x«, s?] . (36) 

Using our parametrization, Zaitsev's boundary condi- 
tions can be solved for the unknown quantities in a 
straightforward way. In the superconducting state we 
define effective reflection and transmission coefficients, 
which we present in Appendix The sum of each gener- 
alized reflection coefficient with its corresponding trans- 
mission coefficient is equal to one. Using these coeffi- 
cients we can write the general boundary conditions for 
the six unknown spin matrix distributions functionson a 
compact form. For the coherence functions we have£3 



rR.A T-}R,A R,A , t~\R,A R.A R,A ryR,A . R.,A r^R,A 

i = ® 7l ' +£>iz ®7 2 =7i ® R ir +7 2 ® D ir > ( 37 ) 

i^R,A r)R,A „ ~ R,A , r\R,A „ -R.A ~ R.A „ ryR.A , ~ R.A „ nR,A fOO\ 

r i = R ii ®7i + D ii ®7 2 ' =7i ®Rir + I2 ® D ir ■ (38) 

Note the intuitively appealing structure of the relations. The outgoing functions are weighted averages of two incoming 

functions. The weights depend on the incoming parameters as well, which reflects the coherence during Andreev 
reflection. The distribution functions have the following boundary conditions 

A^^IZVi* ®A A r , (39) 

At^nVx^ ®A A r . (40) 



JDR 

1 K 



R ii 
K 



Tlxf 8 -in + _M ®Vx* ® —^ 



Rt 
K 
R t 

n 



n R 
u \i 



v 

h R 

V 



V 

V 



Analogous relations, obtained by interchanging the 
subscripts 1 and 2, hold for the other side of the interface. 
The terms proportional to the product VSD — T>(1 — T>), 
are due to particle-hole interference and do not arise in 
the classical limit. Insertion of these equations into Za- 
itsev's boundary conditions shows, that they solve the 
nonlinear problem and eliminate all spurious solutions. 



VI. LINEAR RESPONSE THEORY 

The general linear response theory in terms of the co- 
herence functions and distribution functions was devel- 
oped in Refs. [ p^pi] . Here we give a short review of 
the relevant equations and generalize them for spin de- 
pendent phenomena. For the special case of the diamag- 
netic response see Belzig, Bruder, and Fauchere.EJ We 
assume a small external perturbation and expand g and 
h around the unperturbed solutions. With the replace- 
ments g — > g + Sg and h — ► h + Sh we arrive in linear 
order at the equations 



[e - h, Sg} + ihv f V5g = [Sh, g] ® , (41) 
Sg®g + g®Sg = 0. (42) 

Here the linearized self consistency equations determine 
Sh. For a specially chosen parametrization given at the 
end of Appendices [b| and |c], the linear correction of the 
Green's function, Sg, can be written in terms of the linear 
corrections to the coherence functions, Sy R,A , Sy R,A , and 
linear corrections to the distribution functions, Sx K , Sx K . 

It is convenient to transform from the Keldysh re- 
sponse, Sg K to the anomalous response , Sg a , 

5g a = Sg K ~Sg R ® F eq + F eq ® Sg A , (43) 

with F eq = tanhe/2T. Then, with the definition of the 
anomalous components of the distribution spin matrices, 

Sx a = Sx K +j R ®F eg ®Sy A + Sy R ®F eq ®j A , (44) 
Si a = & K ~ 7 H ® F eq ® Sy A - Sy R ® F eq ® 7 A , (45) 

the spectral response, Sg R,A , and the anomalous response, 
Sg a , are given by 



G 



SfiR,A _ ztOttJ N R ' A Sb I (^^r^+r^®^) (^+ 7 ^®^®7 R ' A ) , ^ frn.A 



(& a — 7 R <g> fi a ®7 A ) — (7 R Sx a — Sx a is 7 A ) 
-(7 B ®& a — <fi a <»7 A ) (<E a - 7 R ® Sx a ® 7 A ) 



Using the anomalous self-energies, fih R A — ( ^ R ' A ^ R ' A \ cfa _ ( <® a ^ a 

~ ^ SA R - A <E K - A y ' ~ ^ -<5A Q -<£ a 

A Q = -Sh R ® F eq + F eq ®Sh A , (48) 



(47) 



(49) 



we define the following short-hand notation for the driv- The spin matrices Sy R ' A , Sy R,A , 6x a , and Sx a are functions 
ing terms in the transport equations, of p/, R, e, t, and satisfy the transport equations 



ihvfV5y R ' A + 2efry r 
ihv/VSj RA - 2eSy T 



(7 R ' A A Jt " 4 + E R - A ) ! 
7 H - A ®<5A i? - A ®7 H - A - 

(^A R ' A + £ B - A ) ( 
1 R - A ®6A R < A ®i R - A - 



= -7 H ®<E a ®7^ 

^ V/ V<S a - ^5 t< S a + (-7 R A n - E R ) 

= -7 J? ®<E a ®7- 4 



, 5y«.^ + fly*.* g (_A«.^«.-4 + 

<E*^ g, 7 «.-4 _ g, ,£ R ' A - 5A R - A 

®& a + & a »(-A A 7 A + S A ) 
+ (SA a ®7 A +7 fl ®(5A a -&] Q , 

»<£ a + <£ a ®(-A A 7 A + S A ) 
+ <5A a ®7 A + 7 R »(5A a -(5E a . 



(50) 
(51) 
(52) 
(53) 



One convenient feature of our parameterization is the 
fact, that the linear response transport equations (|50|)- 
(|53|) decouple for given self-energies. Furthermore, the 
transport equations for 5j R , Sy A , 5x a are stable in direc- 
tion of v/, and the transport equations for Sy R , Sy A , 5x a 
are stable in direction of — v/. This makes a numerical 
treatment much easier than solving the boundary value 
problem for the coupled transport equations (|4l|)-(|42|). 
The It-points for the initial condition correspond to the 
final point or the initial point of the trajectory depending 
on the direction of stability of the transport equation. 



Finally we present the boundary conditions for the co- 
herence functions and for the distribution functions in 
linear response. With an analogous definition of the 
anomalous components of the outgoing distribution spin 
matrices, 



5X° 



SX K + T R 9 F eq (8 St A + ST R (8 F eq 8 f A 



5X a = SX 1 



T R g F eq ® ST J 



5t r ®f, 



eq 



(54) 
(55) 



we obtain the boundary conditions for the corrections 
to the coherence functions and distribution functions, 



<Tf = 


R ll K 




R lr . 


n R 


iVtiy^i 


D R r 






iA R r , 


(56) 


iff = 


R \l x 




R lr 






D R r 








(57) 


SX* = 


R ll 

K 


■giTZSxl c 


n 


V 


%VSx%t 


£>t r 




§ nv&2 i 




(58) 


SX? = 


R ii 
K 




K 


Du 

1 1 

V 




Dt r 




3 TlVtix a 2 i 


5 A A r . 


(59) 
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VII. ANDREEV SPECTROSCOPY AT N-S 
INTERFACES FOR UNCONVENTIONAL 
SUPERCONDUCTORS 

To illustrate the physical content of the introduced 
distribution functions we discuss in this section the An- 
dreev reflection process at an interface between a nor- 
mal metal (subscript '1') and a d-wave superconductor 
(subscript '2'). This ppablem was studied by Blonder, 
Tinkham, and KlapwijkEJ for conventional s-wave super- 
conductors, and was generalized to unconventional super- 
conductors by Bruder.Ei We generalize these calculations 
to include finite impurity scattering and identify features 
which are sensitive to time reversal symmetry breaking 
states. In an Andreev reflection experiment a beam of 
normal quasiparticles with energies, e^, and momenta, 
Pf,b, is injected across the interface into the superconduc- 
tor. Two types of reflections will occur. Part of the beam 
will be regularly reflected at the interface, which amounts 
to a reflection of the quasiparticle's velocity, momentum 
and current, and part will be Andreev reflected. An- 
dreev's retro-reflection is caused by particle-hole conver- 
sion which reverses the velocity but conserves momentum 
and current to very good approximation. Because the 
current is affected quite differently by regular reflection 
and Andreev reflection, a measurement of the current- 
voltage characteristics provides direct information on the 
balance between these two reflection processes. Together 
with a thorough theoretical analysis such measurements 
inform about fundamental properties ot-the supercon- 
ductor such as the symmetry of pairing,c3 the.^sa-p.1 size 
and anisotropics, and interface resonance states.EJcl For 
anisotropic superconductors both the current density in 
the reflected and the Andreev reflected beams will de- 
pend strongly on the direction of the incoming beam, in 
addition to their dependence of the energy of the incom- 
ing quasiparticles. 

The following calculation of Andreev reflection in- 
cludes anisotropic pairing, a finite mean free path in the 
superconductor, a finite transparency of the N-S inter- 
face, the layer of a strongly distorted order parameter 
near the interface, and the effects of the interface on the 
excitation spectrum, in particular the low-energy bound 
states. We consider a layered d-wave superconductor 
with cylindrical Fermi surface and isotropic Fermi veloc- 
ity along the layers. The interface lies perpendicular to 
the layers and we assume, for simplicity, the same Fermi 
velocity in the normal and the superconducting parts of 
the N-S contact. 

The coherence functions 7*, 7J 4 are determined by 
boundary conditions at infinity, Eqs. (C7), JC£|), whereas 
rf , are determined by the interface boundary condi- 
tions, Eqs. (37), pq). For the spin singlet superconduc- 



tor we write 7 



zcr y 7, 7' 



Z<7y7, 



ier„r, and 



r R = i<j y T, where 7, 7, T, T are scalar functions. On 
the normal side the incoming coherence functions 71 , 71 
are zero as a consequence of their zero initial values at 



infinity. Thus, the retarded part of the Green's function 
on the normal side, following from Eqs. j3^), j35|), and 
(U|), has the form 



9i,in 

9l,out 



-ITT 



1 

— 2ia y Ti — 1 

1 2ia y Ti 







-1 



(60) 
(61) 



The nonzero quantities T\ and Ti, describe the proxim- 
ity effect at the N-S interface. The solutions for T% and 
Ti in the normal metal in equilibrium are 



ri(ar,e) = ri(e)e v*e v , 



(62) 
(63) 



where the spatial trajectory coordinate x is measured 
in direction of vj and is zero at the interface, positive 
for Ti and negative for Y\, and n is the lifetime in the 
normal metal. Both amplitudes decay from the inter- 
face towards the normal metal on a scale v/ri. For 
simplicity, we assume in all what follows that the nor- 
mal metal is in the clean limit. The Tomasch oscillation 
factors ,E3 with Tomasch wave length trvf/e, are carried 
by Ti, Ti, whereas 7^ and 72 vary only in the region of 
varying order parameter near the interface and are con- 
stant far away. Similarly, on the superconducting side, 
far away from the interface, the deviations of the out- 
going coherence functions from their homogeneous solu- 
tions, r 2 (x) -T 2 ,hom, r 2 (a:)— r 2 ,/i 0m , carry the Tomasch 



oscillations with wavelength irvf /-\/e 2 — |A| 2 if |e| > |A|. 
In the following all quantities without spatial argument 
refer to their values at the interface. 

In quasiclassical approximation the incoming beam of 
non-equilibrium excitations with energy, tb, and momen- 
tum, p/,6, is described by the "scattering" part of the 
Keldysh Green's function, Ag K = g K — g* q , where the 
equilibrium Keldysh Green's function g^ q is subtracted. 
In the following we assume, for simplicity, a spin un- 
polarizcd incoming beam. The calculations for spin- 
polarized beams pose no new problems bujt are of in- 
terest OBljwior high-field superconductivity,C3 spin-triplet 



pairing Ji-TEj contacts between superconductor^ and mag- 
netic materialSjE-J or spin-active interfaces.tlla The in- 
coming beam is then characterized by unit spin-matrix 
distribution functions Axf and AXf . To obtain a phys- 
ical interpretation of this distribution functions we con- 
sider a solution of Eq. (|l8|) in form of a traveling wave 
with frequency u>, 



Ax*(x,e,t) = Axf (e)e "/ 



(64) 



The corresponding par t of the Keldysh Green's function 
follows from Eq. (|l4| ), and after performing the time 
convolutions, Eq. (Al), we obtain 
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Agf; m (z,e,<) = -2m 



i — (x—vet) 



-f x (e + *)fj(e- 



(65) 





AXf(x,e,i) 



This gives us a very transparent interpretation for the 
processes covered by Axf. The upper left entry describes 
an incoming particle with velocity v/. The lower right 
entry describes an Andreev reflected hole with velocity 
—Vf, coming from the interface due to retro-reflection 
combined with particle hole conversion. The off-diagonal 
components describe particle- and hole-like Tomasch os- 
cillations due to particle-hole coherence. The degree of 
coherence between particles and holes in the incoming 
distribution, Axf, is given by the coherence function Ti. 
This gives a direct physical interpretation of the coher- 
ence functions. Similarly, Ti is the amplitude for An- 
dreev reflected particles due to an incoming hole excita- 
tion beam. On the other hand, the distribution function 
AX* describes an incoherent hole coming from the inter- 
face. This component can be nonzero only if there is an 
incoming hole in the Green's function Agf out or Aggouv 
which we exclude in our scattering boundary condition. 
Thus, the correct boundary conditions for the scattering 
problem take the intuitively appealing form, to allow for 
the incoming particle beam only an incoming distribu- 
tion function, Axf, and for all outgoing channels only 
outgoing distribution functions, AXf , AXf, AXf. All 
other distribution function components are zero. 

In the following we assume a stationary (to = 0) situa- 
tion, where an incoherent beam is injected, which allows 
us to consider the incoming beam spatially homogeneous 
along the trajectory. Furthermore, it is sufficient to solve 
the problem for the distribution function 



Axf = -8ir6e S(e - e b )5(p f - p /ib ) 



(66) 



where p / denotes a unit vector in direction p f , and 5e is 
the energy resolution of the beam. Any other distribution 
of incoming excitations is then given by a linear com- 
bination of such solutions with properly chosen weight 
functions. The current density of the incoming beam is 
jo = 2eNfVfSe. For a current density much smaller than 
the critical current density in the superconductor, one 
can neglect the effect of the beam on the self-consistent 
order parameter and impurity self-energies. 

For the scattering parts of the Keldysh Green's func- 
tion at the normal side we obtain 



Ag£ in = -2th Axf 
Agf out = -2m AX? 



1 



-iOyTi -\Ti\ 

1 




(67) 
(68) 



The vanishing off-diagonal elements of the reflected 
Green's function show that there is no hole admixing in 
the reflected particle beam. 

The boundary conditions for the N-S interface follow 
fromEqs. 



AXf = 72. 
AXf = I? Axf , 

72 



1 + 7272 



1 + 72.7272 



Axf, 



t x =V 



1 + 72-7272 

72 

1 + 72.7272 



AXf = -UV\^ 2 \ 2 Axf 
T 2 = 72-72, 

f 2 = ^72 ■ 



(69) 

(70) 
(71) 

(72) 



The total current densities are given in terms of 
the Keldysh Green's functions via the formula j = 
eNf J (de /8TTi ) Tr(T3Vf Ag K ) B Using the boundary condi- 
tions ©-(72), this gives directly the total current den- 
sities at the interface in terms of the injected current 
density, 



Jl,i 



1+V A 



Jo 

Jl,out 

jo 
jo 

]2,out 

jo 
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= 72. 



1 + 72.7272 

2 



1 + 7272 



nv 



1 + 72.7272 
72| 2 (1 + |72| 2 ) 



V- 



1 



1 + 72-7272 1 

+ I72I 2 



1 + 72.7272 1 2 



(73) 
(74) 
(75) 
(76) 



Here, j\^ n describes the incoming current including the 
excess current, j2,m the regularly reflected current, ji,out 
the regularly transmitted current, and j2,i n describes the 
process where the Andreev reflected holes are regularly 
reflected back to the superconductor at the interface. For 
energies below the gap the transmitted current densities, 
.72, m j j2,out, decay with distance from the interface into 
the superconducting region, where they are converted 
into super-currents. It is straightforward to show the con- 
servation law ji. in +j 2 ,in = jl, out +32, out- EqS. 

hold for general anisotropic and unconventional super- 
conductors, including impurity scattering. The quanti- 
ties 72 and 72 follow from solving numerically their trans- 
port equations, Eqs. (16) and (p7|), with self consistently 
determined self-energies and order parameter. For 








-10.0 -8.0 -6.0 ,-4.0 -2.0 0.0 



-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 



FIG. 2. Order parameter amplitude (left) and local 
density of states at the interface (right) for a interface in 
(100) direction (top) and in (110) direction (bottom). The 
interface is at x — 0, the normal metal extends to x > 0. 
The transmission coefficients for the different curves are 
V =0.1 (full line), V = 0.5 (long dashed), V a = 0.9 
(dashed), and T>o = 1.0 (dotted). The temperature is 
T — 0.3T C , and the mean free path £ = 10£o- 

conventional s-wave superconductors, and assuming a 
step function for the order parameter our formulae agree 
with the results of Blonder, Tinkham and Klapwijk.E] It 
is clear from Eq. (J7j]) that the Andreev reflected beam al- 
ways enhances the current density in the injection beam, 
giving rise to the excess current. The enhancement is 
proportional to D 2 , reflecting the fact that both the in- 
coming particle and the Andreev reflected hole have to 



cross the interface. On the other hand, the current den- 
sity of the conventionally reflected beam, described by 
Eq. (|74|), can be below or above the value 1Z ■ jo- 

The angle resolved density of states at the supercon- 
ducting side of the interface is given by 



N(e,p f ) = N f Re 



1 - 7?-7272 
1 + 7?-7272 ' 



(77) 



The local density of states is given by the Fermi sur- 
face average over this expression. Eq. (77) shows that 
interface bound states are given by the solution of the 
equation 1 + 7?. 7272 = 0. Because the absolute values of 
72 and 72 are in equilibrium always smaller than or equal 
to unity, bound states at an interface can strictly occur 
only for TZ = 1, that means zero transmission. For finite 
transmission the bound states broaden into interface res- 
onances. Impurity scattering further broadens these res- 
onances. In Fig. y we show our self consistent solutions 
for the ci-wave order parameter, A = \^2cos2ip, and for 
the local density of quasiparticle states at the interface. 
For definiteness, we modeled the angular dependence of 
the transmission coefficient for the N-N interface by, 



V{<j>) = 



T>q sin 



K + V a sin z 



(78) 



appropriate for a ^-function potential barrier. Here, <j> 
is the impact angle between incoming trajectory and in- 
terface. The parameters T>q and TZq = 1 — T>q are the 
transmission and reflection coefficients for perpendicular 
impact (<j) = vr/2). The impurity self-energy was calcu- 
lated self consistently in Born approximation with a 
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FIG. 3. Current densities in the injected beam (top panels) and in the regularly reflected beam (bottom panels), 
as a function of the impact angle for 3 energies: et = (full line), 0AA max (dashed) and 1.6A max (dotted). The 
left part of each picture is for a (110) interface, and the right part of each picture for a (100) interface. The left 
picture is for a transmission coefficient T>o ~ 0.5, and the right picture for a transmission coefficient Do ~ 0.9. The 
temperature is T = 0.3T C , and the mean free path I — 10£o- 
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mean free path of I = 10£o- The temperature was chosen 
T = 0.3T C , leading to a maximal gap of A max = 2.29T C . 
For the (100) orientation of the interface the order pa- 
rameter is constant in the superconductor for zero trans- 
mission and is suppressed at the interface for finite trans- 
mission. In contrast for the (110) orientation the order 
parameter is suppressed to zero at the surface foil' = 
and is suppressed to a finite value if T> is nonzero.EU In the 
(100) orientation there is no subgap resonance, whereas a 
zero energy resonance typical for <i-wave pairing at prop* 
erly oriented surfaces is present at (110) orientation.cfl 
Above the maximal gap the density of states is enhanced 
for (100) orientation. There is no such enhancement in 
the density of states at the interface above the gap for 
(110) orientation. 

Figs. H and |] show selected results of our calcula- 
tions of Andreev reflection at a contact between a nor- 
mal metal and a c£-wave superconductor. Our calcula- 
tions are done for T = 0.32^, for three mean free paths, 
£ = 2£o, 10£o, 100£o, and for two orientations of the in- 
terface. Fig. U shows for three energies the dependence 
of the excess current due to Andreev reflection (top pan- 
els) and the regularly reflected current (bottom panels) 
on the impact angle for transmissions T>q — 0.5 (left pic- 
ture) and for transmission T>o — 0.9 (right picture). The 
positions of the gap nodes show up clearly in the Andreev 
reflection amplitude, which breaks down for quasiparti- 
cles transmitted into the nodal directions. The regular 
reflection approaches for the nodal direction the value 



TZ(<f>). The width of this breakdown regions broadens 
with energy. At energies above the maximal gap, A maxi 
the Andreev amplitude approaches zero and the regularly 
reflected amplitude approaches the value The de- 

pendence on the energy of the incoming quasiparticles 
is shown for one representative impact angle in Fig. ^ 
for three values of mean free path, again for transmission 
Vq = 0.5 (left) and for transmission Vq = 0.9 (right). For 
the (100) interface as shown in Fig. [|and |, the behavior 
at low energies is clearly distinct from the behavior for a 
(110) interface. Whereas for a (110) interface the regular 
reflection is suppressed for low energies, it is enhanced 
for a (100) interface. The excess current shows a peak 
at low energies for the (110) interface, but the (100) in- 
terface shows a minimum. The features at the gap edges 
are small for the (110) orientation, but are strong for the 
(100) orientation. And finally, the signal above the gap 
edges is small for a (110) interface but extends well up 
to twice the gap for a (100) interface. 

In the clean limit the zero energy current density of 
the incoming bccim is ji.m /jo = 2 for a (110) inter- 
face, and ii,i„/io = 2(1 + n 2 )/(l +TZ) 2 < 2 for a (100) 
interface; for the regularly reflected beam the zero en- 
ergy limit is ji, ut/jo = for a (110) interface, and 
h,out/jo = 4ft/(l + Tlf > Tl for a (100) interface. The 
values for the (100) interface coincide with the values for 
a conventional isotropic s-wave superconductor, and, are 
in agreement with Blonder et alS3 and Shelankov.E3 




FIG. 4. Current densities in the injected beam (top panels) and in the regularly reflected beam (bottom panels), 
as a function of energy for 3 different mean free path values for the superconductor: £ = 100£o (full line), 10£o 
(dotted) and 2£o (dashed). The impact angle is fyj-K = 0.4. The left part of each picture is for a (110) interface, 
and the right part of each picture for a (100) interface. The left picture is for a transmission coefficient T>o = 0.5, 
and the right picture for a transmission coefficient Do = 0.9. The temperature is T = 0.3T C . The values for the 
maximal gaps at this temperature are A max (£ = 100£o) = 2.13T C (£ = 100£ ), A maw (£ = 10£o) = 12%T c (l = 10£ ), 
A max (£ = 2£ ) = 2.85T C (£ = 2£ ). 
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Explicit values for the zero energy limits at a (100) 
surface are for perpendicular impact ji.m/jo = 1-11, 
ji.out/jo = 0.89 for V = 0.5, and j Mn /j = 1.67, 
ji,out/jo = 0.33 for I? = 0.9. These values agree with our 
numerical calculations for mean free paths £ > 100£o- In 
contrast, the zero energy values for the (110) interface of 
two for the incoming and of zero for the reflected beam 
are very sensitive to impurity scattering. In fact, as can 
be seen from Fig. [|, is the first value reduced to about 
1.2 for half transmission and a realistic mean free path 
of ten coherence lengths, and the second value is larger 
than 0.2 in this case. Also the structures around the gap 
edges for the (100) surface are very sensitive to impu- 
rity scattering. For a mean free path of two coherence 
lengths the Andreev signal is already strongly reduced, 
as our calculations in Fig. |] show. This may explain 
the small signal of only a few percent in many Andreev 
experiments. The different behavior at low energies for 
the regular reflection is the only remaining difference be- 
tween (100) and (110) orientation for mean free paths 
comparable to the coherence length for unconventional 
superconductors. The suppression of the regularly re- 
flected beam at low energies for all angles (except in 
nodal direction), as seen for a (110) interface in the lower 
left panels of Figs. || and ||, is a direct consequence of 
the sign change of the order parameter during reflection 
of quasiparticles. The origin of this effect is the same 
as for the zero energy resonance (and follows from the 
Atiyah-Patodi-Singer theorenfj). Both effects are de- 



stroyed by time reversal symmetry breaking and both ef- 
fects are washed out by impurity scattering. However, in 
contrast to the zero energy resonance, which is not an ex- 
act bound state anymore for finite transmission even for 
zero impurity scattering, the strong suppression at low 
energies of the regularly reflected beam remains a sta- 
ble phenomenon for all transmissions in the clean limit. 
The effect is reduced by finite impurity scattering, and 
in this case it is further reduced if the transmission is 
comparable or smaller than the scattering rate. Thus, 
the zero-energy resonance and the blocking of the reg- 
ular reflection are two complementary phenomena: the 
first one is well established only for interfaces with small 
transmissions, whereas the latter one is well established 
at interfaces where the transmission is not too small. 

The low-energy behavior of the regularly reflected 
beam can be used to prove a sing change of the order 
parameter during reflection of the quasiparticles at an 
interface. Specifically, our results show that for all im- 
pact angles this reflection amplitude is always above the 
normal state reflection, H(4>), whereas for the (110) inter- 
face it is for all directions clearly below IZ(tfi) (the normal 
state reflection can be obtained for a beam with ej, well 
above the maximal gap). 

Finally, we show that the low-energy suppression of the 
regular reflection and enhancement of the excess current 
is a sensitive test for time reversal symmetry breaking 
states. In Fig. [| we show our results for a dominant 
(i-wave coupled to a subdominant s-wave component. 



1.3 i — ■ — i — ■ — i — . — i — ■ — i — ■ — i — ■ — i — ■ — i — . — i — ■ — i 2.0 




e b /A max V* £b /A max V* 

FIG. 5. Current densities in the injected beam (top panels) and in the regularly reflected beam (bottom panels) 
at a (110) interface. Dotted lines are for a temperature T — O.STd, which is above the transition temperature from 
a d to a d + is state. Full lines are for a temperature T = O.lTd, which is below this transition to the spontaneously 
time reversal symmetry broken state. The left picture is for mean free path I — 10£o, and the right picture for mean 
free path I = 100£o- The left part in each picture shows the energy dependence for impact angle 4>/n — 0.348, and 
the right part of each picture shows the dependence on impact angle for e = 0.2A max . The transmission coefficient 
is T>o — 0.2. The subdominant transition temperature is T s ,o = 0.3T C . 
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Below the interface transitioncSO'Ea they couple to the 
spontaneously time reversal symmetry breaking state 
d + isjS'tiZI where the s-wave component is localized in 
a layer of a few coherence lengths near the interface. 
The left picture is for I — 10£o and the right picture 
for £ = 100£o- The coupling strength of the subdomi- 
nant component is characterized by its 'bare' transition 
temperature, T s ^q (in the absence of the dominant com- 
ponent). This transition temperature i& [ted-wced in the 
presence of the dominant component Hj-EZtEa We chose 
T s .o = 0.3Trf, where Td is the transition temperature from 
the normal state to the pure d-wave state. This value of 
T s .o is below the critical .-value for a possible transition 
into a bulk d + is state.c3 Nevertheless, the transition 
into a d + is state localized near the interface is possible. 
According to our calculations the subdominant compo- 
nent is strongly suppressed by finite transmission, so we 
chose a small value, 2?o = 0.2. The dotted curves show 
the current densities of the reflected beams for a tempera- 
ture T = 0.3Td- The system is at this temperature above 
the transition into the time reversal symmetry breaking 
state. Full curves are for T = O.lTd, which corresponds to 
the interface d + is state. As can be seen from Fig. [| the 
suppression of the reflection and the enhancement of An- 
dreev reflection are shifted to negative energies. Due to 
finite impurity scattering, and resulting mixing of differ- 
ent momentum directions, there is also a shadow-feature 
at positive energies. The broadening of the feature itself 
is reduced, leading to a much sharper effect compared to 
the pure d- wave state. The zero energy values for reg- 
ular reflection are changed in the d + is-state to almost 
one. Also shown in Fig. [| is the dependence of the re- 
flected current densities on the impact angle. The small 
dip at (j> — 7r/2 is due to the fact that the energy of the 
incoming particles is above the gap in these directions. 
Below the transition into the d+is state there is a strong 
anisotropy with respect to the interface normal. This ef- 
fect is a consequence of the spontaneous siipercurrents 
at the superconducting side of the interfaced For an in- 
coming particle beam with a projection on the interface 
counter-moving with the current the regular reflection 
is strongly reduced compared to the pure d-wave state, 



whereas the Andreev reflection is enhanced for this case. 
For a beam with a projection co-moving with the super- 
current these effects are absent or inverted. The shift 
of the feature in the energy dependence of the reflected 
current densities is determined by the Doppler shift of 
the quasiparticle spectrum due to the spontaneous su- 
percurrent at the superconducting side of the interface. 
Thus, the sign of the shift for a chosen impact angle can 
be positive or negative dependent on the direction of the 
spontaneous interface currents (similarly the asymmetry 
around the interface normal changes its sign). 

VIII. CONCLUSIONS 

We have presented a new theoretical formulation of non- 
equilibrium superconducting phenomena, including sin- 
glet and triplet pairing, in terms of coherence functions 
and distribution functions. Our central results are equa- 
tions (|l3|)-(|l9|), together with boundary conditions at in- 
terfaces, equations (^)-(^). We used this formulation to 
present the theory for Andreev spectroscopy at interfaces 
between a normal metal and an unconventional super- 
conductor in a transparent way. This formulation allows 
to include disorder in a self-consistent manner. We pro- 
posed an anomalous suppression of the regularly reflected 
quasiparticle beam as a test for time reversal symmetry 
breaking states. This test is especially suitable for not 
too small transmission, where the zero energy interface 
resonances become ill-defined and cannot be used as a 
test for time-reversal symmetry breaking anymore. 
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APPENDIX A: NOTATION 

The noncommutative ® -product is defined in the fol- 
lowing way 



A®B{e 1 t) = e^ 



A(e,t)B(e,t) 



(Al) 



If one of the factors is both independent of e and t, the 
® -product reduces to the usual matrix product. 



For Fourier transformed quantities (t — > lo) we have 

oo 

Auj' du" 



A®B(e,u>) 



2tt 2tt 



■8(u> + to — lo) X 



xi(e+^,w")B(e-^,w') . (A2) 



If A(e, t) = A(e) is independent of t, that means, A is an 
equilibrium quantity, then 
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A®B(e,u)=A(e+^)B{e,w) , 
and, analogously, if B is an equilibrium quantity 

A®B(e,u) =i(e,w)S(e-^) . 
Also we generalize the commutator 

[i,P]® = i®P-P®i. 
The Fermi surface average (■ • -) P ' is defined by 



1 



(2^) 3 |v/(p' 



(A3) 
(A4) 
(A5) 

(A6) 



f> 



were Nf is the total density of states at the Fermi surface 
in the normal state, 



N f = 



(27rn)3|v / (p^)|- 



(A7) 



and Vf(p'f) is the normal state Fermi velocity at the po- 
sition p'y on the Fermi surface 



V /(P/) 



ds(p) 



dp 



(A8) 



p=p> 



Here, e(p) describes the normal state dispersion of the 
quasiparticle band crossing the Fermi level at p'j. 

APPENDIX B: PROJECTORS 

Following Shelankov,0 we introduce the following pro- 
jectors 



(Bl) 



From the normalization condition, g<s>g = — n 2 1, it fol- 
lows that P+ and P_ are projection operators, 



P+ ® P+ = P+ , P_ 8 P- = P- , 

and project orthogonal to each other 



(B2) 



P++P_=l , 
P+®P- = P-®P+ = 0. 



(B3) 
(B4) 



The quasiclassical Green's functions may be expressed in 
terms of P+ or P_ , 



9 = -in (P+ - P-) 
= -in (2P+ - i) = -in (1 - 2P_) 



(B5) 



Equations of motion for the projectors can be extracted 
from the corresponding equations for the quasiclassical 
Green's functions 



[Z-h,P±] + iftv/VP± = 6. 



(B6) 



The Keldysh component of the Green's functions, g K , 
fulfills the relation g R ®g K + g K ®g A = 0. This implies 
P* ® g K ® Pf = 6 and P R ®g K ®P A = 6, leading to 



P R ®g K ®P* + P R ®g K ®P*. 



(B7) 



The value of P R ®g K ®P A does not determine g K 
uniquely. It is possible to add P R ® A + B <g> P A to g K 
with any matrix function A and B without changing 
P+ ®g K ® P_ (similarly for Pf ®g K ® Pf ). One can use 
this property to obtain a proper parameterization of g K 
and to eliminate the unnecessary information in g K . We 
write 



g K =-2ni [P R ®X K ®P A + Pf g F K 



(B8) 



were X K and K K contain only one free spin matrix func- 
tion as parameter. The X K and Y K have to fulfill funda- 
mental symmetry relations, following from the symmetry 
relations for the quasiclassical Green's function, g. From 
the equations of motion of the Keldysh Green's functions 

(ef& - h R ) ®g K -g K ® (ef 3 - h A ) + i%VfVg K = 

-{g R ®h K -h K ®g A ) , (B9) 

we obtain the equations of motion for X K and Y K using 



P R 8 | (er 3 - fr R ) ® A K - X K g (ef 3 - A- 4 ) + ft K + ihv f VX K J ® P 4 = 6 
P* ® ( (era - <g> y K - Y K a (ef 3 - /i A ) - /i' f + ffiv/VF" j ® P A = 6 



(BIO) 
(BH) 
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Tracing these equations properly in the Nambu space 
and respecting the symmetries of X K and Y K , one ob- 
tains two equations for both undetermined 2x2 spin ma- 
trix functions, which parameterize X K and Y K . 

Analogously we proceed for the linear response. From 
the l st -order normalization conditions (J42|) we have 
P R - A 9 5g R ' A <8 P+ A = 6 and P RA ® 5g R - A <g P RA = 6; as a 
consequence the spectral response, Sg R,A , can be written 
as 



6g RA 
=F27ri 



pR.A 



>6X' 



,P* 



Pi 



iP h,a 



(B12) 



Analogously, for the anomalous response the normal- 
ization condition (|4^) leads to Pf ® Sg a ® Pf = 6 and 

P* ® 5g a ® P A = 0, so that 8g a can be written in the fol- 
lowing form, 

Sg a = -2m [Pf ® <it a ® P A + P^ ® 5f a <g> P A 1 . (B13) 



APPENDIX C: PARAMETER 
REPRESENTATIONS OF THE 
QUASICLASSICAL GREEN'S FUNCTIONS 

The projectors Pj* and P R may be parameter- 
ized in the following way by complex spin matrices 
7*(p/, R, e, t) and 7 H (p/, R, e, i) 



P? 



1 



~ R <s (l-7 R ®7 fl ) _1 « (1, 7"), (CI) 



i a (1- 7 R »7")~ ® (l R , I)- (C2) 



Here (1 + argib) is defined by 

(l + a«i6) _1 ® (1 + a®b) = 1 



(C3) 



One immediately proves Pf <g> Pf = Pf , P* ® P* = P* 
and P R ® P* = P* ® P| = 6. A useful identity is 



(l + a®6) ®a = a®(l + &®a) 1 



(C4) 



which may be used to obtain P R + p" = 1 . The unique- 
ness of the projectors is ensured by the symmetry rela- 
tions between the matrix elements of the retarded and 
advanced Green's functions. We may obtain the ad- 
vanced Green's functions either by the fundamental sym- 
metry relation g A = f3(g H yT3 or analogously to the re- 
tarded case using advanced projectors P A = ^(P*)^, 
P A = T 3 (P+) f r 3 



n 



pi 



-7' 
1 



-r 



(1-7 A ®7 A )~ 1 ® (l A , 1), (C5) 
(1 - 7 A ®7 A ) _1 ® ( 1, 7 A )- (C6) 



( 7 K)t 

C6|) into Eq. (|B6|), and using 



holds. Introducing 



Here 7 A = [7 R )^ ; l A 
Eqs. ©, (|c|, ©, 
e <& a + a <& e = 2ea leads to the transport equations for 
ryit.A an( j ^r.a^ gq g ^ ^ ^g|) and (p~7|) , which are generalized 
Riccati differential equations. They are supplemented by 
properly chosen initial conditions. The solutions j R - A 



r 



are introduced into Eqs. (CI), (C2), rtC3), ( 



6Q_to 



obtain the quasiclassical Green's functions, Eqs 
(|l|), viaEq. @. 

The solutions for the coherence functions in a homo- 
geneous singlet superconductor in equilibrium are, 



' rioiri 



i R ' A 

i horn 



signe 



\e R - A \ + \]{e R - A ) 2 + A R - A A R - A 
A R - A signe 



(C7) 
(C8) 



s R,A | + \]{e R - A ) 2 + A R - A A R - A 



where e R ' A = e — (T, R,A — E ff,A )/2. Note that (A R - A A R - A ) 
is proportional to the unit spin matrix and that in the 
clean limit (A R - A A R ' A ) = — |A| 2 . In the presence of a 
constant superflow with momentum p s one has to make 
the replacement e — > e — v/ • p s . 

Using this parameterization the following represen- 
tation for the Keldysh component with hermitian spin 
matrices x K (pf, R, e, t) and x K (pf, R, e, t) is convenient. 
Substituting 



x K 




o o 

x K 



(C9) 



i nto Eq. (|B§| ), using the equation of motion for <? A , Eq. 
(p39|), leads to the transport equations for a:*" and 
Eqs. (18) and (|l9|). Note that e®a — a®e = ihd t a. 
These transport equations have to be supplemented by 
initial conditions. For the Keldysh Green's function Eq. 
@ leads to Eq. Q. 

It is possible to introduce Shelankov's distribution 
functionstZI F and F, which are given by the parame- 
terization 



X' 



F 
F 



F 
F 



(CIO) 



They obey the symmetry relations P(p/,R, e,t) = 
P(- P/ ,R,-e,t)*, F( P/ ,R,e,f) = P( P/ , R, e, i) f - The 
x K and x K are expressed in terms of them in the follow- 
ing way 

x K = (F-7 H «F®fj , 

x K = (P-7 R 0P(57 A ). (Cll) 

Using the introduced parameterization in terms of co- 
herence functions, the transport equation for F has the 
form 
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(ihv f VF + ihdtF)~j R ®(ihv f VF~ihdtF)®f i + 

+ {-Z R ®F + F®i; A + Y, K )- 1 R ®(-£ R ®F + F®S; A - £*)® 7 A + 
— r ) R ® (A R ® F — F ® A A + A K ) + (A R ® F — F ® A A — A K ) ® j A = 



(C12) 



The transport equation for F follows by application of 
the conjugation operation, Eq. (|l2]), to this equation. 
The g K are obtained by introducing Eq. ( Cll ) into Eq. 
(|l4|). The later parameterization is a convenient start- 
ing point for perturbation theory from the equilibrium, 
because in the equilibrium 



eq 



tanh ■ 



2T 



-F, 



eq 



(C13) 



holds and all expression in the braces in Eq. (C12) vanish 
independently. 

Finally we make the connection to our parametrization 
in the linear response. With the choises 



sx R - A 



&y R - J 






LTf R - 



(C14) 



in Eq. flB12j), and 

SX a = 



&c a 




sv a = 





& a 



(C15) 



in Eq. flB13| ), we arrive at equations (^^)-(J53|). With this 
parametrization the linear corrections to the distibution 
spin matrices -y fl ' A , j R,A , x K , x K are given by S~f R,A , dy R,A , 
Sx K , and & K respectively. 



APPENDIX D: REFLECTION AND 
TRANSMISSION COEFFICENTS 

In the superconducting state we define effective reflec- 
tion and transmission coefficients by 



P§ = (1- 7**3 7?), /5g- = (l-7f®7 J R ), (<,j = l,2), (Dl) 

i$ = Kp R 2 g (ftp* + Z>pf 2 ) \ i2* - (ftp* + Pp") ' ® ftp 2 " 2 , (D2) 
ijf,=ftpf 2 ®(ftpf 2 +X>pf 2 ) _ \ £? r = (Kp^+Vp^y 1 ®^, (D3) 

^5 = (7f - 7 2 H ) ® + Z>p£i) ~\ Af r = (ftp 2 R 2 + 2?pf 2 ) ~* ( 7 f - 72 R ), (D4) 

4fi = (7f - 7 2 fl ) ® + -DP2i) ~\ A R r= + Z>pf 2 ) ® (7f - 7?). (D5) 

and Z?f ; = 1 — Rfi, D R r = 1 — R R r D R t = 1 — R R t D R r = 1 — i?f r . Analogously we define these quantities on the other 
side of the interface by interchanging 1 and 2. Advanced quantities are given by the same expressions with the change 
in the superscript R — > A. The following relations are shown to hold 

Al = 1?®^-1?®^, (D6) 
^ = 71^^-72^^ • (D7) 



R 

11 ~ 


R ll 

n 


9 7f- 


V 


3 72 


R 

n — 


R 1i 
n 1 


3 7f" 


f) R 


^72 



As an example we consider the equilibrium spin-singlet r H = ia y T, T R — ia y T, where 7, 7, T, f are scalar func- 
case. In equilibrium the ® -product reduces to a matrix tions. The effective reflection and transmission coeffi- 
product. In this case we can write 7 s = icr y j, 7" = icr y ^, cients are 
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Ri = n 

D X =V 

ki = n 



1 + 7272 

1 + 7^7272 + £>727i 

1 + 727i 
1 + 7^7272 + ^7271 

1 + 7272 

1 + 7^7272 + £>727i 



(D8) 
(D9) 
(D10) 



D X =V 



1 + 7271 
1 + 72-7272 + 7?727i 



(Dll) 



(and analogously for the other side of the interface), 
which fulfill Rj + D 3 ■ = 1 and Rj + Dj = 1 (j = 1,2). 
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